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Abstract 



We study the fluctuation properties of E2 and Ml transition intensities 
among T = 0, 1 states of A = 60 nuclei in the framework of the interacting 
shell model, using a realistic effective interaction for p/-shell nuclei with a 
56 Ni core. It is found that the B(E2) distributions are well described by 
the Gaussian orthogonal ensemble of random matrices (Porter-Thomas dis- 
tribution) independently of the isobaric quantum number T z . However, the 
statistics of the B{M1) transitions is sensitive to T z : T z = 1 nuclei exhibit 
a Porter-Thomas distribution, while a significant deviation from the GOE 

statistics is observed for self-conjugate nuclei (T 2 = 0). 
PACS numbers: 24.60.-k, 24.60.Lz, 21.60. Cs, 21.10.Ky 
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Random matrix theory (RMT) |l] was originally introduced to explain the statistical 
fluctuations of neutron resonances observed in compound nuclei 0. The theory assumes 
that the nuclear Hamiltonian belongs to an ensemble of random matrices that are consistent 
with the fundamental symmetries of the system. In particular, since the nuclear interaction 
preserves time-reversal symmetry, the relevant ensemble is the Gaussian orthogonal ensemble 
(GOE). Bohigas et al || conjectured that RMT describes the statistical fluctuations of a 
quantum system whose associated classical dynamics is chaotic. More recently, a proof for 
this conjecture was proposed by mapping the problem on the supersymmetric non-linear 
sigma model ||. RMT has become a universal tool for analyzing statistical fluctuations in 
chaotic systems 

The chaotic nature of the single-particle dynamics in the nucleus can be studied in the 
framework of the mean-field approximation. The interplay between shell structure and 
fluctuations in the single-particle spectrum can be understood in terms of the classical 
dynamics of the nucleons in the corresponding deformed mean-field potential . However, 
the residual nuclear interaction mixes different mean-field configurations and plays a crucial 
role in determining the nature of the statistical fluctuations of the many-particle spectrum 
and wavefunctions. The statistics of the low-lying collective part of the nuclear spectrum has 
been studied in the framework of the interacting boson model, in which the nuclear fermionic 
space is mapped onto a much smaller space of bosonic degrees of freedom PJlQ|. Because of 
the relatively small number of degrees of freedom in this model, it was also possible to relate 
the statistics to the underlying nature of the mean-field collective dynamics. Since at higher 



excitations additional degrees of freedom (such as broken pairs) become important pT| , it is 
necessary to analyze the effect of interactions on the statistics in much larger model spaces. 
The interacting shell model offers an attractive framework for such studies where realistic 
effective interactions are available (in finite model spaces), and the basis states are labeled 
by exact quantum numbers of angular momentum, isospin and parity ||12|| . 



RMT makes definite predictions for the statistical fluctuations of the eigenfunctions as 
well as the spectrum. The electromagnetic transition intensities in a nucleus are observables 



that are sensitive to the wavef unctions, and the investigation of their statistical distributions 



should complement the more common spectral analysis P,|T0||. B(M1) and B(E2) transitions 
were recently analyzed in 22 Na |L3| in the framework of the shell model and found to 



follow the Porter-Thomas distribution ||14|| , in agreement with RMT and consistent with 



the previous finding of a Gaussian distribution for the eigenvector components \TE\- JTB|. In 
addition, the distributions of E2 and Ml matrix elements did not show any sensitivity to 
spin and isospin fl3|| . 



Most studies of statistical fluctuations in the shell model have been restricted to lighter 
nuclei (A < 40) where complete Ohu calculations are feasible (e.g. sci-shell nuclei). It is of 
interest to investigate how the statistics evolves with increasing mass number. In the present 
paper we study the fluctuation properties of electromagnetic transition intensities in nuclei 
with A ~ 60. We find that the B(E2) distributions are Porter-Thomas, but that the Ml 
statistics is sensitive to isospin. In particular the B(M1) distributions in self-conjugate nuclei 
show a significant deviation from Porter-Thomas. The calculations are performed in the pf 
shell with 56 Ni as a core, i.e. we assume a fully occupied f-j/ 2 orbit and consider all possible 
many-nucleon configurations defined by the 0/ 5 / 2 , 1^3/2 and lpi/2 orbitals. The effective 



interaction is chosen to be the isospin-conserving F5P interaction [Oil. This interaction is 



successful in describing the mass range A ~ 57 — 68. The calculations were performed with 



the shell model program OXBASH [ pO|j . 



Denoting by B{ujL\i — > f) the reduced transition probability from an initial state \i) 
to a final state |/), with u indicating the electric (E) or magnetic (M) character of the 



transition, and 2 L the multipolarity, we have pi 



W-T ITT T T T \ \^T f M ls (uL) - (TjT z 10 1 TfT z ) M iv (uL) 1 2 

BlUjL; Ji±i± z — > Jfl fl z ) = ; ; ; — : . (1) 

v ' } i ' (2Ji + l)(2Ti + l) V ; 

Here Mi S {ujL) and M iv (uJL) are the triply reduced matrix elements for the isoscalar and 
isovector components of the transition operator, respectively. Note that these matrix ele- 
ments depend on Jj, Tj and J/, Tf but not on T z . For AT = transitions (Tj = Tf = T) the 
isospin Clebsch-Gordan coefficient in Eq. (0) is simply given by 



(TT Z 10\TT Z ) = T Z /^T{T + 1) . (2) 

It follows that the isovector component in Eq. ([[]) is absent for self-conjugate (i.e. T z = 
0) nuclei. Thus the statistics of the isoscalar component of an electromagnetic transition 
operator can be inferred directly from T z = nuclei. Consequently, we can test the sensitivity 
of the statistics to the isovector and the isoscalar contributions. 

To study the fluctuation properties of the transition rates, it is necessary to divide out 
any secular variation of the average strength function versus the initial and final energies. 
We calculate an average transition strength at an initial energy E and final energy E' from 
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Ei / B(UL; % -> /) e -( s -^) 2 /27 2 e -(s'-S/) 2 /27 2 
(B(E 1 E')}= '■ — - ^_ {E _ E ^ )2/2 ^^_ (EI _ Ei)2/2i2 . (3) 

For fixed values of the initial (J[,T) and final (JJ,T) spin, isospin and parity, we calculate 
from Eq. (|I|) the intensities B(ZuL;i — > /). All transitions of a given operator (e.g. Ml 
or E2) between the initial and final states of the given spin, parity and isospin classes 
have been included in the statistics. We remark that the energy levels used in (|3|) are the 



unfolded energy levels |23|], characterized by a constant mean spacing. The value of 7 has 
been chosen to be large enough to minimize effects arising from the local fluctuations in 
the transition strength but not so large as to wash away the secular energy variation of the 
average intensity. In the present calculations we used 7 = 4. We renormalized the actual 
intensities by dividing out their smooth part 

B(uL; JjTT z -> J f TT z ) 
Vfi (B(E,E') ' [ ' 

and histogrammed them into bins equally spaced in logy. In RMT there is a large number 
of weak transitions, and we use a logarithmic scale to display small values of y over several 
orders of magnitude. For each transition operator and classes of initial and final states we 
fit to the calculated distribution a x 2 distribution in v degrees of freedom |24| 



P„{y) = (u/2 < y >)^ 2 ^/ 2 - 1 e-^ 2 <^/r(^/2) . (5) 



For v = 1 this distribution reduces to the Porter-Thomas distribution. In systems with 
mixed classical dynamics, it is found that v decreases monotonically from 1 as the system 
makes a transition from chaotic to regular motion |f22f . 

We first examine all E2 and Ml 2+, T -> 2+, T transitions (T = or T = 1) in A = 60 
nuclei. T = states exist only in T z = nuclei, i.e. 60 Zn in our case. However, the T = 1 
states form isobaric multiplets in 60 Co, 60 Zn and 60 Cu. In the latter case we studied the 
statistics in both T z = and T z = 1 nuclei. Because of the vanishing of the isovector Clebsch- 
Gordan coefficient in Eq. (P, the transitions in 60 Zn (T z = 0) are purely isoscalar. For each 
transition operator we sampled 56 2 = 3136 and 66 2 = 4356 intensities for T z = and 
T z = 1, respectively. The calculated distributions (histograms) of the B(E2) (left panels) 
and B(M1) (right panels) 2 + ,T — > 2 + ,T transitions are shown in Fig. [I], and compared 
with the Porter-Thomas distribution (Eq. (||) with v = 1). 

The distributions of transitions within the T = 1 states are shown in the top (T z = 1) 
and middle (T z = 0) panels. Since our interaction is isospin invariant, the energy levels 
of an isobaric multiplet are degenerate, and the spectral statistics of T = 1 states must 
be the same in both T z = and T z — 1 nuclei. We find that the level spacing and A 3 
statistics are in agreement with GOE. However, while the B(E2) distributions are all in 
agreement with the GOE prediction, the Ml distributions show strong sensitivity to T z . 
For T z — 1 the distribution of the Ml intensities is Porter-Thomas, but in self-conjugate 
nuclei we find that the Ml distribution deviates significantly from the GOE limit. The 
dashed line in Fig. |l]d shows a x 2 distribution ([5]) with v = 0.64. Although it does not fit 
well the calculated distribution, the smaller v is consistent with the observed larger number 
of weak transitions as compared with a Porter-Thomas distribution. A similar deviation 
is observed in the Ml distribution for the 2 + ,T = — » 2 + ,T = transitions (see Fig. 
|I]f), while the B(E2) probabilities are distributed according to Porter-Thomas (Fig. |l|e). 
The deviation from the GOE occurs for the Ml transitions in self-conjugate nuclei (e.g. 
60 Zn), where the matrix elements are purely isoscalar. In T z = 1 nuclei both isoscalar and 
isovector components contribute to the Ml transitions. However, since the isoscalar Ml 
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matrix elements are much weaker than the isovector Ml pl |, the latter dominate and the 
distributions are restored to their GOE form. 

We also investigated the statistics of the + ,T = 1 — ► 1 + ,T = 1 transitions and its 
dependence on the number of final states. In the top panels of Fig. ^| we compare the Ml 
distribution (T z = 1) between all 16 + , T = 1 initial states and the lowest 16 1 + , T = 1 final 
states with the corresponding distribution when all 54 final states are taken into account. 
Both seem to give similar distributions, but the statistics is poorer in the first case. 

The Ml strength function is composed of a spin contribution and an orbital contribution, 
each dominating in a different energy region. We find that the orbital part is in close 
agreement with the Porter-Thomas statistics while the dominating spin part shows a slight 
deviation (see Fig. |]c,d). 

Finally, we tested the sensitivity of the statistics to the initial and final spin. In Fig. 
|3] we show the B(E2) distributions for the 2 + — > 4 + and 4 + — > 4 + transitions between 
states with T = 0. The B(E2) distribution for the 2 + — > 4 + transitions is well described by 
Porter-Thomas, as in the 2 + — > 2 + case of Fig. |IJe (in both cases Jj = 2 + ). However, the 
4 + — > 4 + transitions do show a small deviation from the Porter-Thomas limit. This suggests 
that the class of 4 + states is slightly less chaotic. 

In conclusion, we have studied the distributions of B(M1) and of B(E2) transition 
strengths in p/-shell nuclei with A ~ 60. While most of the calculated distributions are in 
close agreement with the Porter- Thomas distribution predicted by the GOE, we find that 
the Ml transitions in self-conjugate nuclei (T z = 0) deviate significantly from the RMT 
limit. For these nuclei the Ml transitions are purely isoscalar and relatively weak compared 
with transitions in T z = 1 nuclei, which are dominated by the isovector component of Ml. 

The authors are grateful for useful discussions with F. Dittes. We thank V. Ponomarev 
for assistance during the preparation of the manuscript. This work was supported in part 
by the Department of Energy grant No. DE-FG-0291-ER-40608. 
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FIG. 1. The B(E2) and £(M1) intensity distributions (histograms) for the 2+,T 2+,T 
transitions in A = 60 nuclei: (a,b) 2+, 1 -> 2+, 1 transitions in 60 Co (T z = 1); (c,d) 2+, 1 2+, 1 
transitions in 60 Zn (T z = 0); (e,f) 2+,0 2+,0 transitions in 60 Zn (T z = 0). The solid lines 
describe the Porter-Thomas distribution (Eq. (||) with v = 1). The dashed lines in d) and f) are 
Eq. (JsT) with = 0.64 and v = 0.54, respectively . 
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FIG. 2. a,b) Distributions of Ml intensities for the + , 1 -> 1 + , 1 transitions in 60 Co (T z = 1) 
with different number of final states included in the statistics (see text). c,d) The distributions of 
the spin and orbital contributions to Ml, for the same transitions as in (b). The solid lines are 
the Porter-Thomas distributions. 
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FIG. 3. The B(E2) intensity distributions for the transitions: a) 2+,T = -> 4+,T = 0; b) 
4+,T = -» 4+,T = in 60 Zn. 
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